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The shape of jamming arches in two-dimensional deposits of granular materials 
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We present experimental results on the shape of arches that block the outlet of a two dimensional 
silo. For a range of outlet sizes, we measure some properties of the arches such as the number of 
particles involved, the span, the aspect ratio, and the angles between mutually stabilizing particles. 
These measurements shed light on the role of frictional tangential forces in arching. In addition, we 
find that arches tend to adopt an aspect ratio (the quotient between height and half the span) close 
to one, suggesting an isotropic load. The comparison of the experimental results with data from 
numerical models of the arches formed in the bulk of a granular column reveals the similarities of 
both, as well as some limitations in the few existing models. 
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I. INTRODUCTION 

Granular materials present interesting and unusual 
physical properties [l|, [2| ■ One of their most salient fea- 
tures is the formation of bridges: stable collective struc- 
tures comprising several grains that can withstand the 
weight above them [3-|a|- Those bridges, or arches, can 
cause the jamming of grains in a fixed configuration that 
is mechanically stable. It has been argued |6| that this 
is a generic property of grain assemblies. In granular 
materials, arches have been proved to be important for 
segregation [7| and non- uniform propagation of forces [8| , 
as well as directly related to the packing fraction and the 
mean coordination number |9l4lll|. In particular, it has 
been found a clear correlation between the steady state 
packing fraction obtained during tapping and the number 
and size of arches in the bulk [12J . 

One most common phenomenon that is a consequence 
of the existence of arches is the arrest of granular flow 
at the outlet of a silo. The jamming of a silo has been 
studied experiment ally both in a two-dimensional con- 
traption (l3rll6| and in a three-dimensional one [17|. In 
both cases, attention was paid to the probability of the 
outpouring being halted by an arch, which depends on 
the ratio between the orifice size and the diameter of the 
grains (if spherical beads are considered). In silo dis- 
charge, for intermediate orifice sizes, arches do not arrest 
the flow; however, they have been shown to be responsi- 
ble for strong fluctuations of the flow rate [18[ . 

In the last years, numerical simulations have been 
widely used to approach this problem. In [19} the effect 
of particle friction and particle size dispersion on arch- 
ing at the outlet was studied. More recently, a simple 
probabilistic model has been shown to capture some of 
the principal aspects of arch formation in hoppers [2Cj . 
In other work [21j , a force analysis of the arches formed 
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following simple rules displays good agreement with the 
jamming probability obtained experimentally Despite 
all these advances in the knowledge of silo jamming, the 
attention payed to the geometry of the arches is scarce. 
To our knowledge, this issue has only been addressed in 
the framework of a study based on the Restricted Ran- 
dom Walk Model (RRWM) to explain arch formation in 
a hopper [15| where some geometrical features were de- 
scribed, and a comparison with the predictions of the 
mentioned model was presented. In addition, there is 
a lack of knowledge about the nature of these blocking 
arches and their precise relationship with the arches that 
are developed within a granular column of deposited par- 
ticles. Intuitively, one could expect that the arches span- 
ning over the orifice are of the same kind as those formed 
in the bulk, but there is no evidence that this is so. 



It is worth mentioning that deciding if a given set of 
touching grains form an arch is a rather complex prob- 
lem. In general, the history of the contacts has to be con- 
sidered J22J. However, the arch that blocks an orifice is 
simpler to study since there are no touching grains below 
the arch. This makes it possible to identify the particles 
that form the blocking arch from a single snapshot. 



In this work, we provide a thorough description of the 
geometry of arches that block the outlet of a silo. To this 
end, we have implemented a 2D setup which allows for 
visual inspection of the arches and a systematic descrip- 
tion of their features. This apparatus is described in the 
following section. We will also compare our results with 
those obtained by numerical simulations of the arches 
within the bulk of a granular column [10|, [22| • This com- 
parison can offer important information about the differ- 
ences and similarities of arches at the outlet and arches 
in the bulk. Some conclusions will be gathered in the last 
section of this article. 




FIG. 1: (a) An arch blocking the outlet of the silo. White 
crosses indicate the centers of the particles that form the 
arch. The beads forming the base (marked with circles) are 
excluded from computations, as explained in the text, (b) 
Example of a bead that hangs from above the equator (the 
angle tj> is larger than 180°). (c) In some cases, the arch forms 
well above the base, with the beads leaning on one or both 
sides of the banks of a quasi-static funnel made of grains, (d) 
Frequency map of the position of the beads for about 20000 
arches that blocked the orifice (units are bead diameters; the 
origin is centered at the exit orifice). Darker areas correspond 
to places where the probability to find a bead belonging to 
an arch is higher; beads at the base have been included. All 
data and photographs in this figure correspond to R = 3.03. 



II. EXPERIMENTAL SETUP AND METHODS 

The experimental apparatus consists of a two dimen- 
sional rectangular silo (500 mm high and 200 mm wide) 
made of two glass sheets between which a thin frame of 
stainless steel is sandwiched. This frame fixes the gap 
between the two glass sheets. At the bottom, there is 
an orifice whose size can be changed at will. The granu- 
lar sample consists of monodisperse stainless steel beads 
with a diameter of 1.00 ± 0.01 mm. In order to give 
an idea of the internal friction of the grains, we have 
measured the avalanche angle, i. e. the slope of the 
plane at which an avalanche develops when a sample is 
slowly tilted from the horizontal. The value of this an- 
gle is 18.5 ± 0.5°, corresponding to a friction coefficient 
of about 0.33 ± 0.02. The grains are poured along the 
whole width of the silo from a hopper at the top. The 
granular deposit within the silo consists of a monolayer 
of particles because the gap between the glass sheets is 
slightly larger than the beads. Finite size effects in the 
lateral dimensions (corresponding to the width of the re- 
cipient) can be neglected as the silo dimensions (width 
and height) are much bigger than the size of the beads 
and the width of the opening at the bottom. A more de- 
tailed description of this device can be found in Ref. [lj| ■ 
The aperture of the exit orifice is given in units of bead 
diameters: R = D/do, where D is the width of the orifice 
and do is the diameter of the spheres. Since do = 1 mm, 
in fact R amounts to the width of the aperture in mm. 

Once the silo is filled, grains pour freely from the exit 
orifice until an arch blocks it. The material flown until a 
jam occurs is called an avalanche. Under proper lighting 
from the back, a still picture (Fig. [T]) of the spheres near 
the orifice is taken (with a standard video camera) when 
the flow is stopped and stored for further analysis. Then, 
the avalanche size is measured — with a scale at the bot- 
tom that collects the grains flown through the orifice — in 
order to validate the measurements with previous results 
(16j |. Finally, the flow is restarted by blowing a jet of 
compressed air aimed at the orifice. The image and data 
acquisition system is computerized and the procedure is 
automated so that the pictures of many blocking arches 
can be registered; the limit is imposed by the capacity 
of the silo, that has to be refilled whenever the level of 
grains falls below a fixed threshold of around 300 mm (1.5 
times the width of the silo) . Seven different outlet sizes 
were explored, and for each one several thousand arches 
were obtained and registered. For some of the orifices, 
about twenty thousand arches were obtained in order to 
increase the resolution in some calculations. 

Even though the particles arrange themselves in a 
monolayer, a small overlap between the spheres in the im- 
ages is possible because the gap between the glass panes 
is necessarily slightly larger than the bead diameter. In 
our device, however, the overlaps detected are smaller 
than about 2 % of the particle diameter. We have per- 
formed several test runs with a granular sample consist- 
ing of round washers, which do not overlap, and we found 



that the appearance of arches is the same. Spheres were 
chosen to allow comparison with previous data that can 
be used as a benchmark (lo |. 

From a digital image such as the one shown in Fig. [TJ 
it is simple to obtain the center of each bead by stan- 
dard techniques of image processing [2J|. In particular, 
we have written a code that involves eroding and dilating 
the image with a disk. Subpixel resolution is attainable, 
but the main limitation for the accuracy in the measure- 
ment of the positions is the overlap of the spheres in the 
images. We can then easily measure the number of par- 
ticles 77 in the arch (those marked with an x in Fig. [TJ 
for instance, r\ = 4 in Fig. QJa)), the span (defined as 
the difference between the horizontal coordinates corre- 
sponding to the centers of the leftmost and the rightmost 
particles marked with white crosses) and the height (the 
vertical distance between the centers of the highest and 
lowest particles marked with white crosses). 



III. ARCH DETECTION 

The definition of arch as a collection of mutually sta- 
bilizing beads requires the knowledge of which particles 
sustain each other in a granular packing. Therefore, the 
protocols designed to identify an arch often involve forces 
between grains or rely on the history of the grains de- 
position [22i . |25| . since, in general, any grain has more 
contacts than needed to make it mechanically stable. 
This information is not easily accessible in an exper- 
iment. This may be one reason for the limited num- 
ber of experimental studies where the arches are studied 
[13l Ha . [25J, whereas there is some literature reporting 
numerical simulations and theoretical studies of arches 
[EL EHHH, [1|. We have found that much simpler con- 
siderations are needed if the focus is shifted to the arches 
that block the exit orifice, leaving all those formed in the 
bulk. 

To identify the blocking arch in a given image from our 
experiment, we take the first line (bottom up) of touching 
beads that span across the orifice (see the beads marked 
in Fig. HJa)). Since our silo has a horizontal bottom, 
beads touching the base are stable per se; therefore, we 
exclude from the arch the two end particles on the base. 
Note that if our silo had the shape of a hopper, the par- 
ticles touching the bottom should be considered as part 
of the arch as they would not be stable per se. In many 
cases, the blocking arch forms up stream instead of at the 
outlet. Beads in the silo may develop a funnel of quasi 
immobile particles that pile up at both sides of the outlet 
and arches may form resting on it [see Fig. [He)]. I n such 
situations the line of touching beads detected will bend 
inwards on one or both sides of the orifice. If the outer- 
most particles do not touch the base of the silo, we define 
the blocking arch as the string of particles that span from 
the leftmost to the rightmost bead in the detected line of 
touching beads, excepting those touching the base of the 
container or forming part of the funnel banks. 




FIG. 2: (Color online)(a) Semilogarithmic plot of the PDF 
of the number of beads r\ that form the arch for different 
values of R, as indicated in the legend, (b) The same results 
as in (a), but the variable is now the number of beads minus 
the diameter of the orifice. The symbols □ and o display the 
results obtained respectively from pseudodynamic simulations 
of hard disks [lOl | and molecular dynamics simulations of soft 
disks with friction 12211. 



IV. NUMBER OF PARTICLES IN THE ARCHES 

One of the most immediate measurements that can be 
obtained from the experiment is the number of beads 77 
that form each arch. In Fig. HJa) the probability den- 
sity function corresponding to different values of R is 
displayed. (The histograms have been normalized so 
that the area below the curve is equal to one). Clearly, 
the number of particles depends strongly on the outlet 
size, because the bigger the outlet the higher the number 
of particles that are necessary to form a blocking arch. 
However, all the distributions display the same trend: an 
exponential decay for high number of particles and a cut- 
off for small number of particles imposed by the outlet 
size. The shape of the number distribution of particles 
in the arches and the exponential nature of the tails im- 



plies that the number of beads per arch has a well defined 
average for a given orifice R. 

It is interesting to take a closer look into the tail of the 
distribution. Based on a pseudodynamic model, it has 
been put forward [lOj that the distribution shows a decay 
sharper than exponential, such as the series represented 
with open squares in Fig.[2]Jb). This has been also found 
in a model of soft disks with static and dynamic friction 
|22j (in these simulations, the dynamic and static friction 
coefficients were set both to 0.5). In probabilistic terms, 
this would mean that the price to pay when a new bead 
is added to an arch would be higher as the arch grows 
larger. In order to assess whether this is the case or 
not, we have obtained the number of beads for around 
20000 arches blocking an orifice of size R = 3.03, which 
corresponds to the data set displayed with filled squares 
in Fig. [2ja). Despite the difficulty of making a definitive 
conclusion, the experimental results seem to indicate that 
the distribution tail is exponential. It should be noted 
that it is impractical to perform more runs because an 
increase of about one order of magnitude in the number of 
arches investigated is needed just to add one data point to 
the plot. With these data, deviation from a straight line 
in a logarithmic scale (meaning an exponential decay) is 
hard to perceive. This would mean that independently 
of the arch size, the probability of adding one bead to an 
arch can be calculated multiplying by some fixed factor, 
as suggested by Mehta et al. [26|. It should be noted 
that in the simulations the distributions depend on the 
packing fraction although in all the cases the decay seems 
sharper than exponential. The assumption made in some 
models, such as those cited in [Tg| , that the decay is 
sharper than exponential, should be reconsidered. 

As mentioned above, from the results shown in Fig. 
Ufa), it is clear that the orifice imposes a cutoff, below 
which no arch can block the exit. Taking as an example 
the results for R = 3.03, arches with -q — 4 seem to be less 
abundant than expected due to this fact. This effect is 
even more accentuated for r\ — 3. It may be thought that 
it is impossible to find an arch of three particles blocking 
an orifice of R = 3.03. However, as there are usually two 
beads touching the silo bottom that are not taken into 
account when computing 77, arches with 77 = 3 are still 
possible for this value of R, although they are scarce. 

We have tried to regroup the PDFs corresponding to 
the different orifice sizes R by rescaling the variable 77 in 
various ways. It finally transpires that just by centering 
all the distribution peaks at the same point regroups all 
the data sets in a single graph, as displayed on Fig. [Hh) ■ 
This in turn implies that there is a linear relationship 
between R and the average of the number of beads: fj = 
1.41 + 1.15 R (see Fig. |3J). It should be noted that a 
linear relationship between the number of beads and R 
is also obtained if the beads forming the base of the arch 
(as explained above) are not discarded. In this case, the 
coefficients of the fit change slightly. We remark here 
that the range of R explored only covers relatively small 
orifices; an extrapolation for larger or smaller orifices of 
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FIG. 3: Average number of beads in the arch that blocks an 
orifice of size R. Error bars correspond to 95 % confidence 
intervals. 



this relationship would be unwarranted. Also note that 
the error bars shown in the plot correspond to confidence 
intervals; the standard deviation is much larger. This 
result (i. e. a linear relationship between rj and R) is 
in good agreement with a previous one that had been 
obtained in a hopper filled with disks for a range of R 
similar to the one explored here [15} . 

The fact that all the normalized histograms, as pre- 
sented in Fig.[2jb), fall on the same curve, indicates that 
the nature of the arches that jam the outlet is the same 
for small and big orifices. In addition, these plots allow 
a straight comparison of the experimental arches with 
those obtained within the bulk through numerical simu- 
lations. In this latter case R — has been taken as there 
is no orifice imposing a cutoff. Despite the difference in 
the tail explained above, the agreement of numerical re- 
sults for the PDF with the experimental ones suggests 
that the arches formed in the orifice (as observed exper- 
imentally) are a representative subset of all the arches 
within the bulk (obtained numerically). Further com- 
ments about this issue are provided in section IVTl 



V. GEOMETRICAL PROPERTIES OF THE 
ARCHES 

Apart from the number of beads that arches comprise, 
there are several features that can be extracted from the 
experimental observations. One of the most interesting is 
the horizontal span [Tj5| . In order to establish a consistent 
comparison with previous numerical simulations, we have 
defined the span as the distance between the abscissae of 
the centers corresponding to the two outermost beads of 
the arch. 

One of the first questions that can be asked when look- 
ing at the arch span is whether the arches formed with 




FIG. 4: [Color online) Normalized PDF (in semilogarithmic 
scale) of the span (s) of the arches formed by rj=5 (magenta), 
6 (green), 7 (blue) and 8 (red) beads obtained with different 
R: A, 2.18; o, 2.69; □, 3.03; V, 3.12; O, 3.61; >, 4.21; <, 4.70. 
Solid lines are the averages over all R for each n. 



a given number of particles display the same span distri- 
bution independently of the size of the outlet that they 
block. Generally speaking, it is found that this is not the 
case (Figs. |U and [S]). But it is true if 77 > R (Fig. HJ), i.e. 
provided that the number of beads is larger than the exit 
orifice, the arches consisting of a given number of parti- 
cles display similar features irrespective of the orifice size 
that they block. In other case, if r\ < R (Fig. [5} , differ- 
ences arise between the span distributions obtained for 
different R, and notably for small spans. This behavior 
can be understood if one considers that, for a given num- 
ber of beads, the smaller spans are strongly affected by 
the cutoff length imposed by the size of the orifice. On 
the contrary, the arches with larger spans are, in all the 
cases, larger than the size of the outlet and hence they 
are not significantly affected by it. 

As the arch span is not severely affected by the orifice 
size (with the condition stated above), it is possible to av- 
erage the data for arches with the same number of parti- 
cles, even if they have been obtained with different values 
of R. We present the averages calculated for arches with 
the same r\ in Fig. 0J Arches in which 77 > R display dis- 
tributions that are close to a Gaussian, whereas smaller 
arches, formed by 2, 3 or 4 particles, display a long tail 
for small span values (Fig. [5]). These results reveal clear 
differences with previous models [T(J, LL3( • Arches of up 
to 4 particles present larger spans than those predicted 
by these models. However, the agreement between mod- 
els and experiment is better for small R (which present 
broader span distributions) since the influence of the ori- 
fice is reduced. Care should be taken when these results 
are compared with those obtained in a hopper [13| , where 
the geometry of the outlet could have a strong influence 
in the arch properties. 

While the span characterizes the size (a direct rela- 



tionship can be established between 77, the span, and the 
height of the arches), the shape can be characterized by 
the aspect ratio: the quotient between half the span and 
the height of the arch. To be consistent with the defi- 
nition of span explained above, the height of an arch is 
defined as the distance between the vertical coordinates 
of the centers corresponding to the highest and the lowest 
beads of the arch. If the aspect ratio is equal to one, then 
half the span is equal to the height. This is obtained, for 
instance, if the arch is semicircular. The aspect ratio is 
large for a flat arch, and small for a pointed arch. 

In Fig. [SI the aspect ratio is represented as a function 
of the number of beads in the arch for different values of 
R. The aspect ratio tends to about one when the number 
of particles in the arch grows. If the beads at the base of 
the arch are included, the asymptotic limit is even closer 
to one. This result supports the hypothesis of semicircu- 
lar arches, as introduced in [15]. Although there may be 
many defects, by which we mean beads departing from 
the angle that would correspond to a perfect semicircle, 
a circular path can be considered a good starting point, 
validating the approach used by the cited authors. Addi- 
tionally, this fact provides a clue about the load sustained 
by the arches. If the arch is optimized to sustain a verti- 
cal load, then it would have an aspect ratio smaller than 
one. On the contrary, an arch optimized to sustain a hor- 
izontal pressure would adopt an aspect ratio bigger than 
one. A semicircular arch (aspect ratio one) is the pre- 
ferred shape to optimize an isotropic pressure [271 l28j . 
Some caution is pertinent here, because a semicircular 
arch could form for any variety of reasons apart from the 
optimization of an isotropic uniform load. Besides, even 
an aspect ratio of one does not guarantee that the arch 
should be strictly semicircular. 

From Fig. [51 we can also observe that big arches have 
approximately the same aspect ratio irrespective of R. 
This indicates that the orifice size does not select a par- 
ticular shape. However, for small arches (small values of 
rj) the plots for different R display different aspect ratios 
which seem to diverge when r\ decreases. Indeed, it is 
clear that the bigger R is, the higher the value of r\ at 
which the divergence takes place. The reason for this be- 
havior is, again, the cutoff imposed by the orifice. When 
the number of particles in an arch is similar to the outlet 
size, the only types of arches that can block the orifice 
are rather flat and hence the aspect ratio is well above 
one. 



VI. ANGLES BETWEEN MUTUALLY 
STABILIZING PARTICLES 

In the previous section, the results of several parame- 
ters that are related with the global shape and size of the 
arches were presented. In order to get further insights on 
the relative position of the particles within the arches we 
have measured the angle <fr between the contacts that sup- 
port each particle in the arch [see Fig. QXb)]. Note that 
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FIG. 5: (Color online) Normalized PDF of the span, s, of the arches formed by (a) 2, (b) 3, and (c) 4 beads obtained with 
different R: A, 2.18; o, 2.69; □, 3.03; V, 3.12; o, 3.61. Dashed lines correspond to the restricted random walk model proposed 
in [13| . Solid lines correspond to arches found in disordered packings of hard disks generated with a pseudodynamic model [10] . 




FIG. 6: (Color online) Aspect ratio A versus the number 
of beads in the arch, obtained for different outlet sizes as 
indicated in the legend. 
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FIG. 7: PDF for the angles associated to the particles form- 
ing an arch, as shown on Fig. UXb). Experimental results 
(circles) and pseudodynamic simulations of hard disks (tri- 
angles) ,10j|. The dashed line indicates (j> — 180°, the angle 
at which a particle is hanging from the equator. Note that in 
the simulations there are no beads beyond cj> = 180°, as they 
are frictionless. 



every particle in the arch, except for the end particles, has 
an associated angle. In Fig. [7] we present the histogram 
of all the angles associated to all the beads forming all 
the measured arches. As we have not found important 
differences between the PDFs when taking into account 
only large (or small) arches, or when considering differ- 
ent spans, heights, aspect ratios, or outlet sizes, all the 
angles for all the beads in all the cases are displayed com- 
bined in the same plot. The probability distribution of 
(j) reveals that there is a large number of particles (about 
17% of all registered beads) with an associated angle 4> 
larger than 180°. These cases correspond to beads hang- 
ing from above the equator, which are stable due to static 
friction since normal contact forces do not contribute to 
the balance of their weight. From these results, we can 
conclude that the effect of static friction is quite relevant 



for arch formation. An example of the important differ- 
ences that arise when tangential friction is neglected is 
observed by looking at the distribution of angles between 
particles obtained from simulations (Fig. [7]). Clearly, the 
beads do not display angles larger than 180°. It must 
be noted that we have observed the same phenomenon 
(a considerable number of particles hanging from above 
the equator) in test runs with washers. This indicates 
that although the gap between the front and rear glass 
panes is a little larger than the bead diameters, the small 
overlapping that can occur is not the cause for the large 
number of particles with <\> > 180°, because these over- 
laps are absent if washers are used. 

An interesting question that arises when considering 
the angles that form the particles in an arch is whether 
there is any kind of correlation among the angles asso- 
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FIG. 8: (a)A probability density map showing the angle cf> 
of the bead i + 1 versus the angle of the preceding bead i 
in the arch, where i is the correlative number of the particle 
in the arch when numbered from left to right. Darker areas 
indicate higher probability; the probability is not normalized, 
and the color bar shows the scale, (b) The PDF for the angle 
4>' = (4>i + 0i+i)/2, which is the average angle between two 
consecutive links. The dashed line indicates <f> = 180°. Note 
the semilogarithmic scale. 



ciated to consecutive particles or not. In Fig. Uta) we 
plot the probability density map of the angle of the bead 
i + 1 versus the angle of the bead i in the arch. It can 
be seen that there is a negative correlation between these 
two quantities. This means that a bead with a large an- 
gle will be likely followed by a bead with a small angle, 
and viceversa. This result implies that the appearance 
of two consecutive beads with large or small associated 
angles is unlikely. Indeed, it seems reasonable that the 
region of </>,; > 180° and 4>i+i > 180° is forbidden as 
two consecutive beads hanging from above the equator 
would conform a very unstable structure. Additionally, 
no angles are detected in the region of </>* < 130° and 
4>i+i < 130° as this configuration would generate an arch 
with a very small span, likely unable to block the outlet. 



The results displayed in Fig. [7] and Fig. [8] (a) seem 
to be at odds with the constraints of the RRWM. In 
particular, the existence of angles <f> > 180 implies that 
the arches are not necessarily locally convex. Hence, the 
condition 9i > 9i+\ (as defined by authors in [15| ) is not 
always fulfilled. Moreover, Fig. [8] hints that the angle 
could be conditioned by the adjacent angles (higher <fii 
means bigger probability of finding a small </>i+i) which 
is against the hypothesis of random angles proposed in 
the RRWM. However, by considering these two contra- 
dictions together, it results that the RRWM is still valid. 
The reason for this lies in the fact that a defect is al- 
ways compensated by the angle associated to its neigh- 
bors. Indeed, if we calculate the mean of two consecutive 
angles <fi' — {<pi + 0,;+i)/2 and we plot its distribution 
(Fig. |U(b)), we can see that the probability of finding 
4>' > 180° is so small that it can be deemed statistically 
insignificant: a clear upper limit in the distribution is dis- 
played at 4>' — 180°. Then, one could relax the condition 
stated for the RRWM to 9i > 9i + 2- Besides, the plateau 
in the distribution around <fr' = 150° — 170° means that 
there is an approximately constant probability for 9, as 
required by the model. Then, although the arch may not 
be convex everywhere due to frictional forces, the RRWM 
is a good approximation as the local concavity is always 
compensated by the neighbors. 



VII. CONCLUSIONS 

In this work we provide a detailed analysis of the geo- 
metrical properties of the arches that block the outlet of 
a 2D silo. The experimental results obtained for different 
outlet sizes reveal that R does not affect the arch shape 
provided that the size of the arch is well above the cut- 
off length imposed by the orifice. We suggest that these 
arches unaffected by the cutoff imposed by the orifice are 
representative of the arches formed in the bulk. Besides, 
we show that, at least for arches with large span when 
compared to the orifice size, an aspect ratio of one (height 
equal to half the span) is preferred. This seems to hint 
that a uniform load (the same from every direction) is at 
play near the orifice. 

The PDF of the number of beads that form an arch 
displays an exponential tail which implies that the prob- 
ability of adding a particle to an arch is constant inde- 
pendently of the arch size [26j. This result contradicts 
previous simulations where it was suggested a sharper 
decay than exponential in 2D [Hj, [22| • 

Additionally, the measurement of the angles between 
consecutive particles in an arch has unveiled the impor- 
tance of static friction forces between grains in the arch 
formation process. In particular, we report that in exper- 
iments there is a noticeable probability of finding angles 
greater than 180°. The existing models of jamming in a 
hopper assume that the only valid angles <f> between the 
particles in the arches are smaller than 180° [13]. When 
friction is taken into account this restriction is fully re- 



leased [13J . However, the correlations displayed by the 
angles of consecutive particles in our experiments may 
alleviate this limitation of the model. As angles bigger 
than 180° are almost ever compensated by neighbors, the 
RRWM would likely be valid if two consecutive angles are 
considered. 

Note that tangential forces are seldom taken into ac- 
count in numerical simulations or models of arch forma- 
tion [3,[a,[2l|. We speculate that the role of static friction 
must be also of great importance in other granular sce- 
narios. In particular, some models of particle deposition 
do no consider this effect 0, [Tfl, [24| and should be revis- 
ited. 
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